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Abstract 

A  geometrical  model  was  developed  to  predict  the  influences  of  solid  grain  size,  pore  size  and  porosity  on  the  triple-phase-boundary  (TPB) 
length  in  electronic  composite  electrodes  of  solid  oxide  fuel  cells.  It  shows  that  the  TPB  length  is  inversely  proportional  to  grain  size  and  can 
be  optimized  by  the  pore  size  and  porosity. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Composite  materials  have  been  investigated  in  the  past 
primarily  because  of  their  enhanced  mechanical  properties. 
For  solid  oxide  fuel  cells,  composite  electrode  materials  are 
of  interest  mainly  because  of  an  increased  density  of  reac¬ 
tion  sites,  thereby  resulting  in  improved  cell  performance. 
The  triple-phase-boundary  (TPB)  in  the  composites  is  the 
active  site  for  electrode  reactions.  Accordingly,  the  length  of 
triple-phase-boundary  plays  an  important  role  in  overall  per¬ 
formance  of  the  cell.  Electrochemical  reactions  in  fuel  cells 
involve  ions,  electrons  and  molecular  species.  In  the  solid  ox¬ 
ide  fuel  cell  (SOFC),  O2  molecules  react  with  electrons  at  the 
cathode  to  form  O2-  ions  which  are  transported  from  the  re¬ 
action  site  by  diffusion  through  the  yttria  stabilized  zirconia 
(YSZ)  electrolyte.  Electrons  are  supplied  via  the  Sr-doped 
LaMn03  (LSM)  electronic  conductor.  The  reaction  occurs  at 
the  triple-phase-boundary  (TPB)  site  where  YSZ  and  LSM 
are  simultaneously  in  contact  with  O2,  i.e.,  an  open  pore.  An 
analogous  situation  exists  at  the  anode  where  Ni,  YSZ  and 
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the  pore  have  a  common  line  of  contact.  It  is  interesting  to  in¬ 
vestigate  the  effect  of  the  ratios  of  phases  and  microstructures 
of  the  composites  on  the  length  of  TPB  and  consequently  on 
cell  performance. 

The  geometry  of  the  anode  has  been  modeled  in  dif¬ 
ferent  ways.  Abel  et  al.  [1]  developed  a  resistor  network 
model  with  Monte  Carlo  simulation.  Results  of  the  calcu¬ 
lation  showed  that  when  an  elementary  reaction  resistance 
Rr  is  not  greater  than  transport  resistance  RQ,  the  maximum 
three-phase-boundary  length  does  not  lead  to  the  largest  ad¬ 
mittance  Y  because  most  of  the  reaction  occurs  at  the  in¬ 
terface.  However,  when  Rr  is  several  orders  of  magnitude 
higher  than  Re,  more  three-phase-boundary  length  is  needed 
to  maximize  the  admittance.  Tanner  et  al.’s  model  [2]  as¬ 
sumed  that  an  electrolyte  material  in  an  electrode  may  be 
viewed  as  “tree-like”  and  protruding  from  the  dense  elec¬ 
trolyte  surface  with  electrocatalyst  particles  spread  over  the 
surface  in  a  connected  network.  His  results  suggested  a  sig¬ 
nificant  benefit  from  making  the  microstructure  (including 
porosity)  as  fine  as  possible.  However,  a  contradictory  result 
was  derived  from  Maggio’s  model  [3]  which  indicated  an  op¬ 
timum  range  of  pore  radii  necessary  to  reach  the  maximum 
limiting  current  density.  Smaller  pore  sizes  may  result  in  poor 
limiting  current  density. 
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In  this  paper,  relationships  between  the  length  of  the  TPB 
versus  grain  size,  pore  size  and  porosity  are  modeled  based 
on  geometric  analysis,  providing  some  insight  into  the  con¬ 
tradiction  in  the  literature. 


2.  Modeling 

This  model  assumes  an  electronic  composite  which  con¬ 
tains  two-solid-phases  and  a  gas  phase  (pore),  the  solid  grains 
being  spherical  and  packed  as  shown  in  Fig.  1.  The  elec¬ 
tronic  and  ionic  conducting  phases  are  present  in  equal  vol¬ 
ume  amounts. 

Two  grains  of  two- solid-phases  with  an  identical  diame¬ 
ter  Ds  form  a  neck  contact  of  diameter  aDs  (a  <  1)  between 
them.  The  exposed  surface  area  of  the  solid  grains  is  equal  to 
the  surface  area  of  the  pore  phase.  The  number  of  two- solid- 
phase-boundaries  (A^2pb)  can  be  obtained  from  the  following 
equation  when  the  solid  phases  have  equal  volumes  and  dis¬ 
tribute  uniformly: 

AU°  -  N2pbAx  =  Ap  (1) 


and 

=  V4  =  Vs  =  614 

s  Vs°  (4/3)tt(Ds/2)3  jzDI' 

N' i  = _ Y* _ 

P  vp°  (4/3)7T(Z)p/2)3  tz  D3’ 

A°s=nDl  Ap  =  ApTVp  =  ^jr~, 

uv 

Vs  +  Vp=l,  Ax  =  (|)  D3(  1  -  v/l  -  a2) 

where  A^pb  is  the  number  of  two-phase-boundary  per  unit 
volume;  Ns  the  number  of  solid  grains  per  unit  volume;  Nv 
the  number  of  pores  per  unit  volume;  A^  the  surface  area  of 
one  solid  grain;  Ax  the  surface  area  disappeared  due  to  the 
contact  between  two-solid-grains;  Ap  the  total  surface  area 
of  the  pore  phase  per  unit  volume;  Ds  the  diameter  of  a  solid 
grain  =  grain  size;  Dv  the  diameter  of  a  pore  =  pore  size;  Vs° 
the  volume  of  one  grain;  Vs  the  volume  of  solid  per  unit 
volume;  Vp  the  volume  of  pores  per  unit  volume  =  porosity. 


Fig.  1.  Equal-sized  grains  of  two-solid-phases  in  contact  (one  an  ionic  con¬ 
ductor  and  the  other  an  electronic  conductor)  form  a  two-phase-boundary 
with  diameter  aDs  (a  <  1). 


Fig.  2.  3D  graph  of  the  relationship  between  TPB  length  vs.  grain  size  and 
porosity  when  the  porosity  is  30%  (a)  and  the  pore  size  is  5  |xm  (b). 


The  length  of  the  three-phase-boundary  is  given  by  the 
following  equation: 


^3pb  2<.  N2pb(7C0tDs)Ap  — 


72 q[Dp  -  (Dp  +  ps)Vv]Vv 

£>2 £>2(1  _  ^/l  _a2) 


Fig.  3.  3D  graph  of  the  relationship  between  negative  log  of  the  TPB  length 
vs.  pore  size  and  porosity  when  the  solid  grain  size  is  5  |xm. 
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Fig.  4.  Relationship  between  the  TPB  length  and  grain  size  of  solid  phase  when  porosity  and  a  are  set  to  30%  and  0.7,  respectively.  The  pore  size  varies  as 
shown. 


3.  Results  and  discussion 

Eq.  (2)  represents  the  triple-phase-boundary  length  as  a 
function  of  Ds,  Dp  and  Vp.  The  TPB  length  is  inversely  pro¬ 
portional  to  the  solid  grain  size  Ds,  and  the  relationship  be¬ 
tween  the  TPB  length  and  Dp  and  Vp  is  schematically  shown 
in  Figs.  2  and  3  when  a  is  approximated  to  0.7.  There  is 
little  variation  of  TPB  length  over  realistic  value  of  a ,  and 
therefore  the  effect  of  a  is  not  considered  further. 


The  graphs  in  Fig.  2  clearly  show  that  smaller  grain  size 
greatly  increases  the  TPB  length  no  matter  how  the  pore  size 
and  porosity  change.  However,  the  TPB  length  decreases  sig¬ 
nificantly  and  TPB  domain  shrinks  when  the  grain  size  in¬ 
creases.  For  larger  grain  size,  the  TPB  only  covers  the  domain 
of  larger  pore  size  and  lower  porosity. 

Fig.  3  indicates  that  a  smaller  pore  size  results  in  high  TPB 
length  only  when  porosity  is  low  and  further  reducing  the  pore 
size  would  decrease  the  length.  Overall  TPB  length  decreases 


Fig.  5.  Relationship  between  the  TPB  length  and  grain  size  as  a  function  of  pore  volume  (porosity).  The  pore  size  is  5  [xm. 
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Fig.  6.  Relationship  between  the  TPB  length  and  pore  size.  The  porosity  is  30%. 


with  increasing  pore  size  but  larger  pores  can  tolerate  much 
wider  porosity  range. 

The  2D  relationship  between  the  TPB  length  and  solid 
grain  size,  pore  size  or  porosity  will  be  discussed  in  more 
detail  below. 

First  we  examine  the  relationship  between  the  TPB  length 
and  grain  size  which  is  plotted  in  Fig.  4  with  porosity  fixed 
at  30%  and  a  at  0.7. 

Overall,  the  TPB  length  dependence  exhibits  three  regions 
(Fig.  4).  The  triple-phase-boundary  length  decreases  expo¬ 


nentially  with  increasing  grain  size  when  the  grain  size  is 
small,  for  example,  the  TPB  length  drops  by  two  orders  of 
magnitude  when  the  grain  size  changes  from  0.5  to  2  [xm. 
Hence,  controlling  the  solid  grain  size  is  extremely  impor¬ 
tant  in  maximizing  the  TPB  length.  With  increasing  grain 
size,  the  TPB  length  exhibits  a  relatively  lower  decrease. 
Larger  pore  sizes  extend  this  stage  of  behavior.  In  the  third 
stage,  the  TPB  length  decreases  rapidly  with  increasing  grain 
size.  A  worst  case  scenario  occurs  when  pores  are  much 
smaller  than  the  grains  and  the  TPB  length  drops  precipi- 


Fig.  7.  Relationship  between  the  TPB  length  and  pore  size  at  different  porosities  of  0.2,  0.3  and  0.5  and  solid  grain  size  of  5  [xm. 


X.  Deng,  A.  Petrie  /  Journal  of  Power  Sources  140  (2005)  297-303 


301 


Fig.  8.  Relationship  between  optimal  pore  size  and  porosity  for  different  solid  grain  sizes  of  1  and  10  pan. 


tously.  Large  pores  can  tolerate  a  wider  distribution  of  grain 
sizes  before  the  TPB  length  suffers  a  severe  drop.  Gener¬ 
ally  speaking,  a  finer  microstructure  results  in  a  longer  TPB, 
which  is  in  agreement  with  Tanner’s  model  [2].  A  nanos- 
tructured  composite  would  have  good  potential  for  SOFC 
application. 

Fig.  5  shows  the  relationship  between  TPB  length  and 
grain  size  at  different  porosities  when  the  pore  size  is  5  pun. 
The  TPB  length  again  decreases  with  increasing  solid  grain 
size.  Note  that  at  a  high  porosity,  the  TPB  length  would  be  less 


tolerant  of  large  solid  grains.  In  other  words,  a  higher  porosity 
cannot  support  a  large  TPB  length  as  grain  size  increases. 

The  pore  size  dependence  of  the  TPB  length  is  shown  in 
Figs.  6  and  7.  It  indicates  that  TPB  length  increases  sharply 
with  increase  of  pore  size  and  then  reaches  a  maximum  be¬ 
fore  decreasing.  The  decrease  is  at  a  much  lower  rate  than 
the  increase,  e.g.  when  the  grain  size  is  5  |jim,  the  TPB  length 
increases  by  orders  of  magnitude  until  the  pore  size  reaches 
3  [xm  but  remains  near  the  same  value  after  3  |xm.  Conse¬ 
quently  there  is  an  optimal  pore  size  to  maximize  the  TPB 


Fig.  9.  Relationship  between  the  LTPB  and  porosity  at  different  pore  sizes  of  1,  5,  10,  and  20  |xm.  The  solid  grain  size  is  5  pm. 
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Fig.  10.  Relationship  between  the  TPB  length  and  porosity  at  different  solid  grain  sizes  of  1,  5,  10,  and  20  |xm  and  a  pore  size  of  5  pan. 


length,  which  is  in  agreement  with  the  results  of  Maggio’s 
model  [3]  and  results  [4].  The  optimum  pore  size  (Z)ppt)  at 
which  the  TPB  length  is  maximized  shifts  to  a  higher  value 
as  grain  size  and  porosity  increase.  Dppt  can  be  obtained  by 
differentiating  Eq.  (2): 

d(^3pb)  q 

dDp 

from  which  we  obtain: 


£>°pt  =  2  Ds 


Yl 

1  - 


From  Eq.  (3),  the  Z)pPt  is  directly  proportional  to  the  solid 
grain  size  Ds  and  porosity  Vp  when  Vp  is  very  small,  but 
approaches  infinity  as  Vp  goes  to  1  (Fig.  8).  In  other  words,  a 
larger  pore  size  would  optimize  the  TPB  length  as  the  grain 
size  and  porosity  increase. 

In  the  SOFC  case,  a  target  of  30%  porosity  is  usually  main¬ 
tained  in  the  electrodes.  If  Vp  =  30%  is  substituted  into  Eq.  (3), 

we  obtain  Dppt  =  (6/7 )DS.  It  is  suggested  from  Figs.  6  and  7 
that  it  is  safe  to  maintain  the  pore  size  larger  than  6/7  of  grain 
size  since  the  TPB  length  drops  significantly  if  the  pores  are 
smaller  than  Z)ppt. 


Fig.  11.  Relationship  between  optimal  porosity  and  pore  size  at  different  solid  grain  sizes  of  1  and  10  pm. 
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Finally  Figs.  9  and  10  show  the  dependence  of  the  TPB 
length  on  porosity.  The  TPB  length  is  also  optimized  by  con¬ 
trolling  the  porosity  at  different  pore  or  grain  sizes.  Previous 
work  on  SOFC  performance  measurement  [4]  supported  the 
effect  of  porosity  on  the  TPB  length.  The  optimal  porosity 
VpPt  at  which  the  TPB  length  is  maximized  can  be  derived 
from  Eq.  (1)  in  the  same  way  as  above: 


2(T)p  +  T)s) 


It  is  indicated  from  Eq.  (4)  that  when  the  pore  size  is  much 
smaller  than  solid  grain  size,  VpPt  increases  linearly  with  pore 
size  but  reaches  a  limit  when  the  pore  size  becomes  much 
larger  than  the  solid  grain  size  (Fig.  11). 

The  curves  of  TPB  length  versus  porosity  are  symmetric 
about  FpPt,  which  implies  that  deviation  from  VpPt  in  either 
direction  would  decrease  the  TPB  length  equally. 


4.  Conclusion 

As  an  important  parameter  of  a  three-phase  (two- solid 
and  one-pore-phases)  electronic  composite,  the  three-phase- 
boundary  length  is  a  function  of  size  and  volume  fraction  of 
the  three  phases.  A  finer  solid  grain  size  increase  the  TPB 
length  exponentially,  while  the  pore  sizes  and  porosity  have 
optimal  values  at  which  the  TPB  length  can  be  maximized 
and  the  optimal  values  depend  on  grain  and  pore  size  and 
porosity.  The  curves  of  TPB  length  and  pore  size  are  asym¬ 
metric.  The  TPB  length  drops  significantly  faster  on  one  side 
of  the  maximum  than  the  other  side.  However,  the  depen¬ 
dence  of  the  TPB  length  on  porosity  is  symmetrical  about 
the  optimal  porosity.  Any  deviation  of  porosity  from  the  op¬ 
timal  decreases  the  TPB  length  equally. 

We  can  combine  the  results  of  this  analysis  with  the  con¬ 
straints  on  the  fabrication  of  solid  oxide  fuel  cell  electrodes. 
The  electrodes  are  sintered  from  micron-sized  powder.  The 
resulting  grain  size  is  approximately  5  jxm  but  could  vary 
from  1  to  10  |xm  or  more.  Pore  size  is  non-uniform  as  some 
pores  originate  from  inter-particles  space  and  some  from 
graphite  which  is  added  as  a  pore  forming  agent.  We  gen¬ 
erally  expect  a  pore  size  in  the  range  of  1-10  [xm  but  larger 
pore  size  could  be  produced  if  desired.  The  target  porosity 
is  30%  with  an  allowable  range  of  20-40%.  There  is  limited 
open  porosity  below  20%  and  poor  crush  strength  above  40% 
porosity. 

From  our  geometrical  analysis  aimed  at  maximizing  the 
TPB  length,  we  conclude  that 


1.  A  pore  size  of  5-10  jxin  diameter  is  ideal  when  the  grain 
size  is  1-10  |xm  (at  30%  porosity).  A  1  [xm  grain  size  has  a 
TPB  length  2-3  orders  of  magnitude  greater  than  a  10  pm 
grain  size. 

2.  A  fine  grain  size  of  1-2  pm  and  high  porosity  (30%  or 
greater)  offer  the  best  combination  for  maximizing  TPB 
length.  On  the  other  hand,  for  a  grain  size  distribution  of 
3-10  pm,  less  than  30%  porosity  is  optimum. 

3 .  The  TPB  length  has  a  low  dependence  on  pore  size  beyond 
a  critical  limit.  For  30%  porosity,  the  critical  ratio  of  pore 
to  grain  size  is  6/7.  Very  poor  performance  is  predicted 
below  the  critical  ratio. 

4.  The  TPB  length  is  not  highly  sensitive  to  changing  poros¬ 
ity  near  the  optimal  value.  For  example,  for  a  5  mm  grain 
size,  the  TPB  length  changes  by  less  than  a  factor  of  2 
between  10  and  40%  porosities. 

5.  An  ideal  combination  of  parameters  is  predicted  to  be  a 
pore  size  of  1-2  pm,  a  grain  size  of  1-2  pm  and  a  poros¬ 
ity  of  30-60%  in  order  to  maximize  the  number  of  TPB 
sites.  From  practical  consideration,  it  would  be  difficult 
to  limit  grain  growth  to  such  small  size  during  sintering 
and  operation  and  provide  adequate  gas  diffusion  through 
such  fine  pores. 
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